The four height variables of the Abehan sandpile model 
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We study the four height variables in the Abehan sandpile model. We argue that the four variables 
are not represented by the same operator along closed boundaries, or in the bulk. Along open bound- 
aries, we calculate all n-point correlations, and find that there, all height variables are represented 
by the same operator. We introduce dissipative defect points, and show that along open boundaries 
they are represented by the same operator as the height variables. We show that dissipative defect 
points along closed boundaries, or in the bulk, have no effect on weakly allowed cluster variables. 

PACS numbers: 05.65.+b,45.70.-n 
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The Abelian sandpile model (ASM) proposed by Bak, 
Tang and Wiesenfeld produces power laws without any 
fine-tuning of parameters, and thus potentially provides 
an explanation for how power laws can arise in nature fl| . 
Since its introduction, the ASM has been used to analyze 
a diverse range of systems — see for a review. 

The ASM is an extraordinarily simple mathematical 
model (see for a description). However, significant as- 
pects of the ASM are still unknown. While the height 
one variable is well understood, the higher height vari- 
ables (two, three and four) are not. We approach this 
problem in this paper primarily by looking at correlation 
functions along or near open and closed boundaries. We 
also investigate the role of dissipative defect sites. 

The correlation functions are computed with an ele- 
gant method, introduced by Majumdar and Dhar, who 
used it to calculate unit height probabilities and correla- 
tions 3] . The recurrent states of the ASM can be mapped 
to spanning trees drawn on the sandpile lattice ^ |^ . A 
spanning tree is a connected, directed, acyclic graph, such 
that every site has a path leading to the root, which is 
"off the edge of the sandpile" (i.e. the "site" connected 
to all open boundaries). The site i is said to be a pre- 
decessor of the site j if the path from i to the root goes 
through J. The probability for a site i to have height h is 
equal to the probability that in a spanning tree, exactly 
h — \ oi its nearest neighbors will be predecessors of i 
(NNP's, for nearest-neighbor predecessors) [^. 

Certain height configrations in the ASM correspond 
to local restrictions on the spanning tree. For example, 
the condition for a site to have height one is equivalent 
to a spanning tree condition that it be disconnected from 
three of its neighbors. Any probabilities corresponding to 
local restrictions on the spanning tree (known as weakly 
allowed cluster variables) can be calculated as finite- 
dimensional matrix determinants with the Majumdar- 
Dhar method ^ |^. These calculations require use 
of the lattice Green function. 

Spanning trees are described by the central charge —2 
logarithmic conformal field theory (c = —2 LCFT), and 
so the question arises as to how variables in the ASM 
should be identified with operators in the LCFT. Mahieu 



and Ruelle, using the Majumdar-Dhar methods, calcu- 
lated correlation functions for a number of height config- 
urations, and showed that they agree with LCFT correla- 
tions, with appropriate field identifications 7]. However, 
these methods do not allow the calculation of probabil- 
ities or correlations for bare higher height variables, be- 
cause they are associated with predecessor relationships, 
which are nonlocal; a site i can be a predecessor of a 
neighbor j by a long path which goes far from either i or 
j. Priczzhev has developed O-giwph. methods to find the 
bulk probabilities for higher height variables 6]. How- 
ever, the bulk correlations for higher height variables re- 
main unknown. 

The sandpile boundaries can either be open, where 
grains of sand can fall off the edge during topplings, or 
closed, where they cannot. Ivashkevich has shown that 
along either boundary type, the nonlocal arrow diagrams 
associated with higher height probabilities and two-point 
correlations can be written as linear combinations of lo- 
cal arrow diagrams jS] . He found that all boundary two- 
point correlation functions, between all height variables, 
fall off as and thus argued that all height variables 

correspond to the same LCFT field. We reanalyze his 
calculations and results below. Mahieu and Ruelle pre- 
sented other evidence that the higher height variables 
are identical to the unit height variable in the scaling 
limit; however, they also pointed out that this identifi- 
cation appears inconsistent with LCFT operator prod- 
uct expansions (OPEs) J^. We argue here that analysis 
of boundary correlation functions, and heights far from 
the boundary, indicate that the height variables should 
not all receive the same field identification along closed 
boundaries, or in the bulk. However, along open bound- 
aries, they are identical in all n-point correlations; addi- 
tionally, dissipative defect points along open boundaries 
receive the same field identification as the height vari- 
ables. 

Closed: We define, for all n-point correlation functions 
along closed boundaries, 
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where "c" stands for "closed," the a^'s are heights, is 
the height at position x along the boundary, and Pai,c is 
the constant probability for a site on the closed boundary 
to have height a^. Note that along closed boundaries, 
sites cannot have height four. Ivashkevich found 
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While these correlations all have the same power law, the 
coefficient in /c(l, 1) is negative, while the coefficients of 
/c(2, 2) and /c(3, 3) are positive. These sign differences 
indicate that the three height variables are in fact not 
all represented by the same field operator. Furthermore, 
the coefficients of /c(a, a') in do not factorize (into 
KaKa'), as would be expected if all three height variables 
were the same up to rescaling. 

As a check, we have rederived all results of . (There 
is a misprint in the result for /c(3,3) in j^]-) While we 
agree with the results of @, the analysis there appears 
to have several errors. For the two-point correlations 
between and rj, it divides ASM states into sets Sab, 
where the states of Sab are allowed with heights hr- > a 
and hr > b, but forbidden otherwise. However, not all 
ASM states fall into such sets. For example, there are 
states allowed when (hn, hr^) = (1, 2), and allowed when 
{hr-,hr-) — (2, 1), yet forbidden when {hr-, hr-) — (1, 1). 
However, the relationship in between the Sat and the 
spanning trees is also not quite correct, and these errors 
end up largely, but not entirely, cancelling. The spanning 
tree representation is somewhat surprising. The natural 
assumption (made in [^) is that all spanning trees where 
ri and rj each have one NNP contribute to the 2-2 corre- 
lation in the ASM. However there are some such spanning 
trees that do not; for example, trees where and rj each 
have one NNP, and neither is a predecessor of the other, 
but a neighbor of is a predecessor of r^, and a neigh- 
bor of rj is a predecessor of , contribute not to the 2-2 
correlation, but to the 2-3 correlation. Such cases can be 
written as linear combinations of closed loop diagrams, 
and then calculated using a generalized Kirchoff theo- 
rem 0. Luckily, these graphs fall off as l/(xi — 
leaving the results of ^ unaffected. 

We have calculated all closed boundary three-point 
correlation functions that include at least one unit height 
variable. Some results are 
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The correlation functions must satisfy identities result- 
ing from the fact that the height probabilities at any site 
must sum to one. For example, X]a3=i fc{ai, 0,2, 03) = 0. 
As a check, we have verified that all other three-point 
correlations that include at least one unit height vari- 
able satisfy all such identities. All calculated two- and 
three-point correlations are consistent with bulk correla- 
tions identifying -(2(37r - 8)/tt'^) [dOde) with the height 
one variable, (6(7r - A)/tt^) {dOde) +{1/27^)6 with the 
height two variable, and (8/7r2) {8080)- {I / 271)98^9 with 
the height three variable. There are no 9's, since fields 
near a boundary of a (L)CFT can be thought of as 
purely holomorphic; CFT boundary correlations behave 
like bulk correlations, where the antiholoniorphic parts of 
the field act like holomorphic pieces at mirror locations 
across the boundary 0, ^| . It is also consistent to 
make the substitution 9^9,9^ —9 in all field identi- 
fications above, as the c — —2 CFT is symmetric under 
this transformation. {9 and 9 are variables in the c = —2 
LCFT. See for a brief description of the c = -2 LCFT, 
and references therein for a more complete treatment.) 

We have not been able to calculate all terms contribut- 
ing to /c(2,2,2), but note that if complications arising 
from graphs such as those (to be) discussed below Figs.Q] 
and [2 are ignored, we obtain 



/c(2,2,2) 



(24 - 57r)(-576 + 3847r - Glir^) 
ATr^{xi ~ X2Y{xi - X3)^{x2 - x^y 



+ 



(8) 

which is consistent with the field identification above for 
the height two variable above. 

We have also introduced a dissipative defect site on the 
closed boundary. At a normal site of a closed boundary, 
a site topples whenever it has more than three grains, 
losing three grains, and sending one grain to each of its 
three neighbors. On an open boundary, sites topple when 
they have more than four grains, and thus send an extra 
grain off the edge. We add a "dissipative defect point" 
on the closed boundary, where the site can now have up 
to 3 -I- fc grains, and where k grains fall off the edge with 
each toppling. Previously, dissipation has been added 
uniforml y in the bulk, where it breaks criticality |l2l| . 
(But see [l3j for dissipation along a line.) We have shown 
that if the defect-free closed boundary Green function is 
Gc. 0(^1 J ^2), then a defect at d changes the Green function 
(dropping constant terms) to 



G,{rl.,r2) = Gc,o(rl, rs) - G,.o{rl,d) - Gc,o(d, ri) (9) 
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The modified Green function is independent of the value 
of k. (This is reasonable, since regardless of the amount 
of dissipation, the defect provides the only possible span- 
ning tree route to the root.) Using this new Green func- 
tion, for a defect at a; = on the boundary, we find 



/c(l) = 
/c(2) - - 
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The height one probability is unaffected, to all orders^ at 
all points, in the bulk and the boundary, by the defect 
(except at the defect itself). J2l=i /c(«) = 0, so /c(3) = 
—fc{2). Eqs. 1)1011 l|l have been numerically confirmed. 

These results indicate that a dissipative defect on a 
closed boundary is represented by a dimension opera- 
tor. (Note that Mahieu and Ruelle have identified uni- 
form dissipation in the bulk of the ASM with the integral 
of a dimension operator 0.) The defect has correla- 
tions with the height two and three operators, but not 
with the height one operator, again pointing to differ- 
ent field identifications for the different height variables. 
Since they differ along closed boundaries, they must also 
differ in the bulk, since CFT boundary operators are de- 
rived from OPE's of bulk operators 

Finally, we have found the correlation function of n 
unit height variables, at positions Xi {i = l,2,...,n) 
along the boundary. This requires local arrow constraints 
at 3n vertices of the ASM, and thus the calculation of a 
3n-dimensional matrix determinant. The matrix is di- 
vided into 3 by 3 block submatrices, such that the diag- 
onal blocks are all identical, and the off-diagonal blocks 
all have the same form. A rotation makes the matrix 
diagonal in 2 out of every 3 rows (and columns). The 
universal part of the correlation function is found to be 



— ] det (M) , 



where M is the n-dimensional matrix 
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This is the same as the n-point correlation of 
-(2(37r - 8)/7r2) [dBdO) in the bulk, confirming the unit 
height identification below Eq. (j?)). In 0, the unit height 
variable in the bulk was associated with dOdO dOdO. 
Open: The correlation functions are simpler along open 
boundaries. We define the operators 
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Pa, op is the constant probability for a site along an open 
boundary to have height a (already found in Q), and the 
K„ are normalization factors: 
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Using methods similar to those used to find Eq. H12|) . we 
find that the n-point open boundary correlation function. 



(</)ai (a;i)0a2 (2:2) ■• • 0a„ (a^n)) 
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is equal to det(M). (These results reproduce and extend 
the one- and two-point functions found in |^.) Eq. I|16|) 
is independent of the a^'s, showing that the four height 
variables, upon rescaling, do all receive the same field 
assignment {0989) along open boundaries. Apparently 
the different height variables correspond to different fields 
in the bulk, and remain different along closed boundaries, 
but become the same along open boundaries. 

We again add dissipative defects. At the site d = {x, y), 
we increase the toppling condition by /c > 0, so that k 
grains of sand are dissipated with each toppling. We 
assume that all open defect sites are y — 0(1) from the 
boundary (at y = I). Then the new Green function is 

k ^ ^ 

-Y^iGopA^I, d)Gop,oid,r2) (17) 

1 + fcGop,o(a, a) 

Gopflifi ,f2) is the defect- free open boundary Green func- 
tion. Eqs. H17I) and lO are different because the Green 
function between nearby points is 0{1) on an open 
boundary, but OihiL) on a closed boundary, where L 
is the distance to the nearest open boundary. We have 
generalized Eq. (|17|l for multiple open dissipative defects. 

We define (f)^ {d; k) as the operator corresponding to the 
addition of a defect of strength k at d = (cc, y), and then 
the multiplication of all correlation functions containing 
(f>5{d; k) by 7r(l-|-fcGop,o('^: d))/{ky'^). Then the connected 
n-point correlation function in Eq. (|16|l is still given by 
the connected terms of det(M), even if some of the at are 
now equal to 5, and regardless of the (different) values 
of the fc's at the various defects. So the addition of a 
local dissipative defect near an open boundary is, like 
the height variables, represented by d9d9. 

The open boundary is much more tractable than the 
closed boundary for several reasons. For calculating one- 
point functions on any boundary, when we write the non- 
local arrow diagrams as linear combinations of local ar- 
row diagrams, we use the equivalence of certain nonlocal 
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FIG. 1: Two graphs identical only in one-point functions 

Boundary Boundary 




FIG. 2: A long path 



arrow diagrams. For example, Ivashkevich pointed out 
that Figs. ^ and ^ are equally likely because we can 
switch from^ to^ by reversing arrows along the long 
path 0. However, Figs.^ and^^ are not equivalent if 
they are embedded in a correlation function, which has 
arrow conditions in distant parts of the lattice, because 
reversing the long path can change these distant arrows. 
This produces extra terms, such as the graph in Fig. El 

The graph in Fig.Elcan be written as a linear combina- 
tion of closed loop graphs. It, and all analogous graphs, 
turn out to fall off faster than l/(xi — 0:2)* for the open 
case, and thus make no contribution to any universal 
correlation functions. On the other hand, along closed 
boundaries, the graph in Fig. and other analogous 
graphs, fall off as i/{xi — ^2)^, complicating matters (al- 
though, in the end, the l/{xi— parts of these closed 
loop paths end up cancelling). 

Also, the Green function decays as along open 

boundaries (unlike the closed and bulk cases, where it 
grows as ln(a;) 0|). This allows trace formulae found 
by Mahieu and Ruelle for 2- and 3-point correlations 
to be generalized for open boundary n-point functions, 
simplifying matters. 

Bulk: Analysis of the higher height probabilities well in 
the bulk, at distances y ^ 1 from a boundary provides a 
further argument that the higher height variables should 
not all receive the same bulk field identification. Suppose 
that they do all have the same field identification, up to 
rescaling. Then we would know both that 

Paiy) = Pa,B + -| + . . . (18) 

y 

for all a = 1 . . . 4 (based on results for a = 1 in [l3|). 
and also that all height a-height a correlations in the 
bulk would have to have the same sign (negative 0]). 
Upon rescaling the height variables to give them (neg- 



ative) unit norm, the Ca are rescaled to Ca- Based on 
general (L)CFT arguments, the Cj, should be universal 
numbers, independent of a 

m O El- (See in for 
an explicit demonstration of this, in the context of the 
ASM.) Then the original (unrescaled) Co must all have 
the same sign. However, the four Ca's cannot all have the 
same sign, since Y^'^=iPa{y) — 1 for all y. So, contrary 
to our assumption, the four height variables must have 
different bulk field identifications. 

A dissipative defect in the bulk modifies the Green 
function exactly as in Eq. , and again does not modify 
unit height probabilities at any sites other than the de- 
fect. However, modest numerical simulations show that 
higher height probabilites are changed near the defect. 
This parallels the the closed case, again indicating differ- 
ent field identifications for the higher height variables. 

More generally, we can show that no height configura- 
tion probabilities that can be calculated by the removal 
of a set of bonds in the ASM (i.e. the weakly allowed 
cluster variables, and their correlations) are affected by a 
dissipative defect in the bulk or closed cases. This lack of 
correlations suggests that the analysis of weakly allowed 
cluster variables in 0, while impressive, should not be 
expected to give a fully representative picture of the field 
structure of the ASM. 

After the bulk of this work was completed, we were 
informed of independent, then unpublished, calculations 
of Eqs. (I5| 7 |l ■ in the massive case, by G. Piroux and P. 
Ruelle 16]. Discussions with G. Piroux and P. Ruelle 
then led us to the field identifications below Eq. Q, and 
they further corrected an error we had made in these 
same equations. This work was supported by Southern 
Illinois University Edwards ville. We thank V. Gurarie 
and E. V. Ivashkevich for useful discussions. 
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